5 PAGES, 2 FIGURES, ACCEPTED FOR PUBLICATION IN APJL 
Preprint typeset using I4TgX style emulateapj v. 10/09/06 



O 

o 

(N 

1—5 



o 



> 

o 
o 

On 
O 



LONG-TERM STABILITY OF NON-THERMAL X-RAY MODULATION 
IN THE GAMMA-RAY BINARY LS 5039 

Tetsuichi Kishishita'-^, Takaaki Tanaka\ Yasunobu Uchiyama^-*, and Tadayuki Takahashi'-^ 

5 pages, 2 figures, accepted for publication in ApJL 

ABSTRACT 

We report on long-term stability of X-ray modulation apparently synchronized with an orbital period of 3.9 
days in the 7-ray binary LS 5039. Recent observations with the Suzaku satellite in the year 2007, which covered 
continuously more than one orbital period, have provided us with detailed characterization of X-ray flux and 
spectral shape as a function of orbital phase. Motivated by the results from Suzaku, we have re-analyzed the 
X-ray data obtained with AS'CA, XMM-Newton, and Chandra between 1999 and 2005, to investigate long-term 
behavior of LS 5039 in the X-ray band. We found that the modulation curves in 1999-2007 are surprisingly 
stable. Even fine structures in the light curves such as spikes and dips are found to be quite similar from one 
orbit to another The spectral characteristics observed in the past are consistent with those seen with Suzaku 
for some orbital phase segments. We suggest that magneto-hydrodynamical collisions between the relativistic 
outflow from a compact object and the stellar wind from the O star explain the clock-like non-thermal X-ray 
emission over eight years through remarkably stable production of high-energy particles near the binary system. 
Subject headings: X-rays: individual (LS 5039) — X-rays: binaries — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Gamma-ray binaries, a new sub-category of X-ray binaries 
that radiate very high energy ( VHE) 7-rays, have been attract- 
ing much attention since such sources were detected by mod- 
ern air Cherenkov telescopes. Up to now, Te V 7-rays are de- 
tected from four binary systems; LS 5039 (lAharonian et al.l 
I2005ah. PSR B1259- 63 ( lAharonian et "ai] |2005bl). Cvgn us X- 
1 (lAlbert et ai1l2006h . and LS I h-61 303 (lAlbert et alll2007l 
lAcciari et al.ll2008h . While a compact star powering the sys- 
tem is known to be a young radio pulsar in PSR B 1259-63, a 
black hole in Cygnus X-1 (still an evidence of TeV 7-ray de- 
tection), it remains controversial in the case of LS 5039 and 
LS I H-61 303. Also, HESS J0632-t -057 is pointed ou t as a 
candidate to new gamma-ray binary dHinton et al.ll2009l) . 

The LS 5039 system was designated as a high mass X-ray 
binary based on the ROSAT sW-sky survey data by lMotch et al.l 
(II997I) . The stellar companion is a bright star with V=11.2 
and its spectral type is 06.5V((f)). Th e distance to the bi - 
nary is estimated as c/ = 2.5 ±0.1 kpc (ICasares et al.ll2005h . 
The compact object is orbiting around the stellar compan- 
ion with an orbital period Poi-b = 3.9 0603 ± 0.00017 day s 
and with eccentricity e = 0.35 ±0.04 (ICasares et al.ll2005l) . 
Milliarcsecond-scale collimated outflows ext ending up to 
~1000 AU were detected by radio observations (iParedes et al.l 
|20S[2002; Ribo et al. 2008). 

Recent observations by H.E.S.S. revealed that LS 5039 is a 
source of Very High Energy (VHE) 7-rays a nd hence, is able 
to acce lerate particles to multi-TeV energies (lAharonian et al.l 
l2005al) . Moreover, H.E.S.S. detected an or bital modulation 
of the VHE emission with its orbital period (lAharonian et al.l 
l2006h . The Hght curve as a function of orbital phase has a 
peak around inferior conjunction (when the compact object 
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appears between the optical star and the observer). 

Most recently, new X-ray observations with Suzaku have 
enabled us to fully characterize the flux and spectral shape 
as a function of phas e over one and a half orbital periods 
dTakahashi et al.ll2009l) . The Suzaku observations showed that 
the X-ray spectrum at each orbital phase can be well fit by 
a power law and that photon indices vary within F = 1 .45- 
1.60 depending on orbital phase. Based on the Suzaku and 
H.E.S.S. data, it was concluded that the X-ray emission is 
likely to be synchrotron radiation from relativistic electrons 
which are also responsible for the VHE 7-rav emission via 
inverse Compton scattering dTakahashi et al. 2009). 

Now that the X-ray characteristics of LS 5039 for more than 
one orbital period has been uncovered by Suzaku, it is crucial 
to compare the temporal and spectral behavior as a function 
of orbital phase between the Suzaku data and past X-ray data. 
Whether or not the orbital modulation is stable during the past 
decade should provide us with clues to understand the emis- 
sion mechanism and the nature of the compact object in the 
LS 5039 system. In this Letter, we report on the discovery 
of the long-term stability of the X-ray orbital modulation, by 
combining the archival ASCA, XMM-Newton, and Chandra 
data with the recent Suzaku data. 

2. DATA REDUCTION 

LS 5039 has been observed several times in X-rays for 
limited period with ASCA , XMM-Newton, and Chandra (see 
iBosch-Ramon et al] 12007). The log of those observations as 
well as the Suzaku observation is given in Table[T] Data reduc- 
tio n and analysis method s for the Suzaku data are described 
in iTakahashi et al.l d2009l) . Barycentric correction was per- 
formed for all the data in the following analysis. 

2.1. ASCA 

The scientific payload of ASCA consisted of two co- 
aligned detector systems, the Gas Imaging Spectrometer 
(GIS; Ohashi et al. 1996) and th e Solid-state Imaging Spec- 
trometer (SIS: lBurke etaOl991l) . The GIS and SIS detectors 
covered the energy range of 0.7-10 keV and 0.5-10 keV, re- 
spectively. Considering the low-background and large effec- 
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TABLE 1 
Results of Observations 



Date 




Mission 


Observation ID 


MJD 


Exposure (ks) 


Phase 




pa 


Hux(l-10 keV) ' 


X W) 


1999-10-04 


ASCA 


47001000 


51455.5 


24.4 


0.30-0.49 


6.4 ± 1.7 


1.55 ±0.12 


9.59 ± 0.28 


5.7(13) 














0.30-0.40 


6.4 (fix) 


1.55 ±0.10 


8.73 ± 0.43 
















0.40-0.49 


6.4 (fix) 


1.54 ±0.07 


10.42 ± 0.36 




2003-03- 


-08 


XMM-Newton 


0151160201 


52706.3 


10.3 


0.52-0.55 


7.2 ± 0.6 


1.51 ±0.07 


9.20 ±0.17 


8.2 (6) 


2003-03- 


-27 


XMM-Newton 


0151160301 


52725.8 


10.3 


0.53-0.56 


7.0 ± 0.5 


1.45 ±0.07 


8.93 ±0.17 


5.9 (6) 


2004-07- 


-09 


Chandra 


4600 


53195.1 


11.0 


0.67-0.71 


7.9 ±0.8 


1.49 ±0.06 


12.33 ±0.28 


5.2 (6) 


2004-07- 


-11 


Chandra 


5341 


53197.5 


17.9 


0.28-0.34 


7.9 ±0.8 


1.53 ± 0.05 


7.33 ±0.17 


4.5 (10) 


2005-09- 


-22 


XMM-Newton 


0202950201 


53635.8 


15.1 


0.46-0.51 


6.8 ± 0.4 


1.43 ± 0.05 


10.36 ±0.15 


8.9(12) 


2005-09- 


-24 


XMM-Newton 


0202950301 


53637.8 


9.8 


0.99-0.02 


7.0 ± 0.9 


1.57 ± 0.10 


5.99 ±0.14 


3.4 (9) 


2007-09- 


-15 


Suzakii 


402015010 


54352.7 


203.2 


0.00-1.47 


7.7 ±0.2 


1.51 ± 0.02 


8.07 ± 0.06 





Quoted errors are at the 90 % confidence level.'' Nh units are 10 
is calculated from the modeled orbital light curve. 

live area at high energies, we used only the data from the GIS 
detectors. During the observation the GIS was operated in a 
normal PH mode, and standard screening procedures were ap- 
plied for the analysis. The source photons were accumulated 
from a circular region with a radius of 4.'5 while background 
region was chosen in the same field of view with the same ra- 
dius and an offset of 18' from the source. The averaged source 
and background count rates in the energy range of 1-10 keV 
are 0.12 cnt s~' and 0.02 cnt s~', respectively, for each GIS 
detector. In order to improve statistics, the source photons of 
both GIS detectors (GIS2 and 3) were co-added. 

2.2. XMM-Newton 

For the XMM-Newton observations, we analyzed data from 
the European Photon Imaging Ca mera (EPIC), which consists 
of two MOS (iTurner et al.ll2001h and one PN ( IStriider et all 
IMlCCD arrays. All the observations were performed with 
the Medium filters; the observing modes were Prime Small 
Window (time resolution ~6 ms) and Prime Partial W2, for 
the PN and MOS cameras, respectively. Reduction and anal- 
ysis of the data were performed following the standard proce- 
dure using the SAS 7.0.0 software package. After discarding 
high-background time intervals, the source photons were ac- 
cumulated from a circular region with a radius of 40". Since 
the fluxes measured by the EPIC PN instrument show sys- 
tematically ~ 10-15% smaller values in comparison with the 
fluxes measured by the EPIC MOS, we used only MOS data 
and co-added source photons of both detectors (MOS 1 and 2). 
The averaged source and background count rates in 1-10 keV 
are 0.45 cnt s~' and 0.02 cnt s"', respectively, for each MOS. 

2.3. Chandra 

We analyzed data from Advanced CCD imaging Spectrom- 
eter (ACIS; Garmire et al. (2004])) for Chandra observations. 
Because of the superb angular resolution of the X-ray mir- 
ror (~ 0."5 for on-axis observations), the Chandra data could 
be easily affected by pile-up, particularly when a source is 
observed on-axis. The pile-up effect causes underestimate 
of source flux and derivation of harder spectral index. In 
four data sets in the Chandra Data Archive, LS 5039 was 
observed on-axis for two of them (ObsID: 2730, 6259), and 
was observed off-axis for th e other two (ObsID: 4600, 5 341). 
In the previous analysis by iBosch-Ramon et al] (|2007|) . the 
spectra of the on-axis observations actually showed harder 
photon indices (F ~ 1 ■ 1 ) c ompared with the Suzaku results. 
iBosch-Ramon et al] (l2007h noted that their results are likely 
to be affected by the pile-up effects. We also analyzed the 
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Fig. 1. — Orbital variations of photon index and unabsorbed flux in the 
energy range of 1-10 keV. Each color indicates the data of XMM-Newton 
(blue, cyan, and green), ASCA (red), and Chandra (magenta). Black filled 
and open circles represent the Suzaku observation. Fitting parameters are 
shown in Table [T] 

same data sets and concluded that the pile-up effects are sig- 
nificant and difficult to be accurately restored. We therefore 
made use of the data from the two off-axis observations to 
avoid large systematic errors. In the data analysis, the Chan- 
dra Interactive Analysis of Observations software package 
(CIAO 3.4) was used to extract the spectrum, following the 
standard procedure given in the Chandra analysis thread. The 
source photons were accumulated from a circular region with 
a radius of 45" while the background region was chosen in the 
same field of view with the same radius and an offset of 2'.5 
from the source. The averaged source and background count 
rates in 1-10 keV are 0.40 cnt s"' and 0.02 cnt s"', respec- 
tively. 

3. ANALYSIS AND RESULTS 
3.1. Spectral Analysis 

iTakahashi et alJ (l2009h reported that the Suzaku spectra of 
LS 5039 are well fit with power laws and the photon indices 
vary within F = 1.45-1.60 depending on orbital phase (see 
Figure [T]i. We performed spectral fitting for the all archival 
data to compare with the Suzaku results. 

All the spectra are well represented by absorbed 
power-law functions (using the wabs model in XSPEC, 
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Fig. 2. — (a) Orbital light curves in the energy range of 1-10 keV. (Top) Suzaku XIS data with a time bin of 2 ks. Overlaid in the range of 4> = 0.0 to 2.0 is the 
same light curve but shifted by one orbital period (open circles). (Bottom) Comparison with the past observations. Each color coiTesponds to XMM-Newton (blue, 
cyan with each bin of 1 ks, and green with each bin of 2 ks), ASCA (red with each bin of 5 ks), and Chandra (magenta with each bin of 2 ks). Fluxes con'espond 
to unabsorbed values. The blue solid lines show periastron and apastron phase and red dashed hues show superior conjunction and inferior conjunction of the 
compact object, (b) Closeup in 1.2 < < 1.8. 



iMorrison & McCamnionI Il983h . In Table [T] we show the 
resulting fitting parameters, which agree with the literature 
(jBosch-Ramon et al. 2007; Martocchia et al. 2005). The pho- 
ton indices obtained for the ASCA, XMM -Newton, and Chan- 
dra are within F = 1.45-1.55, and consistent with the Suzaku 
values within statistical errors. In addition, the A^h values ob- 
tained are generally consistent with the Suzaku best-fit values 
although the statistical errors are relatively large. 

In Figure [T] we plot the photon indices and flux in the 
range 1-10 keV as a function of orbital phase. Here, the 
orbital phase is calculated with the period of 3.90603 days 
and = 0.0 with reference epoch 7b (H JD - 2400000.5 = 
51942.59) taken fromCasares et al.i (l2005h . Typical cross cal- 
ibration uncertainty for photon index is .1, estimated from 
independent fits of Cr ab spectrum F ~ 2. 1 (iKirsch et al.l2005t 
iKokubun et al.ll2007l) . We fixed the photoelectric absorption 
column density at A^h = 0.77 x 10^^ cm"^ for phase-resolved 
Suzaku spectra, which is obtained from the phase-averaged 
Suzaku spectrum. Also, phase-resolved ASCA spectra were 
fit using A^H frozen at the phase-averaged value. As seen in 
Figure [T] the photon indices obtained in the past observations 
follow the tendency seen in the Suzaku data, where the in- 
dices become smaller (F ~ 1 .45) around apastron and larger 
(F ~ 1 .60) around periastron. Te V 7-ray e mission presents 
similar trends in flux and photon indices dAharonian et al.i 
I2OO6I) as those of X -rays. Remarkably, the X-ray flux shows 
almost identical phase dependency between the Suzaku and 
the other data sets as seen in the bottom panel of Figure[T] The 
similar te ndency in fluxes and photo n indices were also pre- 
sented in iBosch-Ramon et al.l (l2005l) . although the data were 
background contaminated and this prevented a proper photon 
index determination. The results above indicate that the over- 
all orbital modulation has been quite stable over the past eight 
years. Variability on shorter timescales is studied in the fol- 
lowing subsection. 

3.2. Temporal Analysis 

The Suzaku light curves bv lTakahashi et al.l (l2009l) revealed 
variability on short time scales of Ac^ ^0.1. To investigate the 



long-term behavior of the X-ray modulation, we compared the 
light curve by Suzaku and those obtained in the past observa- 
tions. To directly compare the light curves obtained with the 
different detector systems, we need to convert detected count- 
ing rate to absolute energy flux for each bin of the light curves. 
For all the data other than the Suzaku, we assumed power-law 
spectra with parameters shown in Table [T] and converted the 
counting rate of each time bin to power-law flux (unabsorbed 
flux). We adopted variable bin widths to have equalized errors 
for different data sets. For the Suzaku data, since the spec- 
tral parameters, particularly photon indices, are significantly 
changed during the observation, we converted the observed 
counting rate assuming the power-law parameters obtained 
The left panel of Figure |2] shows the resulting light curves 
in the energy range of 1-10 keV. Also, a magnified plot of the 
light curve is shown in the right panel of Figure|2] The phase- 
folded X-ray fight curves from the ASCA, XMM-Newton, and 
Chandra observations are in remarkable agreement with those 
from Suzaku. More surprisingly, short-time variability seen in 
the Suzaku 2007 data can be recognized in the previous 1999- 
2005 data. For instance, a small peak around (p = 0.70 is ev- 
ident both in the Suzaku and Chandra data although the two 
observations are performed almost three years apart. When 
the Suzaku data points in 0.6 < (p < 0.72 excluding the peak 
structure are approximated to a linear function, the deviation 
of this peak structure from the continuum corresponds to 5.6 
fj significance. The flux drop around (p = 0.51 is also consis- 
tent between the Suzaku and XMM-Newton data within sta- 
tistical errors. The deviation of the drop structure from an 
extrapolated linear function determined from the Suzaku data 
points in 0.4 < < 0.5 is 7.4 <t significance. The ASCA data, 
which were obtained 8 years before the Suzaku observation, 
also reproduce the Suzaku light curve not only for the over- 
all flux increase but also for the small peaks around = 0.40 
and (f> = 0.48. Since the X-ray spike around 4> = 0.70 is located 
around inferior conjunction (when the compact object appears 
between the stellar companion and the observer), it seems to 
be due to a rather geometrical than physical effect in the sense 
that it is connected to a direction to the observer In this re- 
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gard, the Doppler effect is a possible way of boosting or mod- 
ul ating the X-ray flux alon g the orbital phase, as pointed out 
in iBogovalov et all (l2008h . It should be noted that there are 
slight differences in the phase-folded light curves. Although 
a small spike at = 0.70 in the Suzaku light curve agrees well 
with the Chandra light curve as seen in the left panel of Fig- 
ure 12] it is evident in the right panel that the Chandra spike 
lasts longer than the Suzaku one by ~ 0.01. Also, only 
the second orbit data of Suzaku shows a dip-like structure at 
(/) = 0.35, while there are no such structures in the other data 
including the first orbit data of Suzaku. 

The largest systematic errors in flux of each time bin are 
due to uncertainties in the cross calibration between each de- 
tector system. Although cross calibration for a source with 
hard spectrum like LS 5039 has not been performed exten- 
sively, the discrepancy in the absolute flux is thought to be 
less than 20% in the current calibration uncertainties (e.g. see 
http:www.iachec.orgiachec_2008_meeting.html for the cur- 
rent cross calibration issues). This number is rather large 
if compared to the amplitude of the small structures seen in 
the light curves. However, the agreement of the similarity 
of shapes between Suzaku and other observatories would not 
simply be a chance coincidence in most cases since the agree- 
ment is seen for more than one observation by Chandra and 
XMM-Newton. In order to quantitatively estimate the simi- 
larity of shapes, we interpolated the Suzaku light curve and 
compared with each orbital segment's light curve. The inter- 
polated values are calculated with the Suzaku light curve of 
5 ks time bins and a spline function, and values are derived 
from the difference between the interpolated and each light 
curves. Each value is pesented in Table[T] The model light 
curve is acceptable in the confidence level of 10-90% from 
test for each segment. We note that the issue whether the 
source is really periodic or not at X-rays, despite the appear- 
ance, is still not solved statistically in periodicity analysis. 

4. DISCUSSION 

We have shown that the orbital modulation of LS 5039 is 
surprisingly stable over the past eight years. It is remarkable 
that not only the overall modulation but also the fine structures 
in the light curve agree well between the data from Suzaku and 
those from other X-ray observatories. 

Considering the long-term stability, it is quite difficult to 
attribute the X-ray emission of LS 5039 to the emission from 
hot plasma around an accretion disk, which is often observed 
from black hole binaries in low/hard state, since such emis- 
sion is generally expected to show unstable an d unpredictable 
time variability dLewin & van derKlisll2006b . Although the 
jet formation and the production mechanism of non-thermal 



radiation in the jet are still unclear, it might be difficult to 
explain the emission from LS 5039 as synchrotron emission 
from a jet, since jet activities are thought to be connected 
to accretion activities and disk instabihti es affect the jet for- 
mation and radiation ('Fe nder et alj|2004l) . Among currently 
proposed models for LS 5039, a scenario in which the ra- 
diation is related to an ultrarelativistic pu lsar wind (iDubusI 
120061: ISierpowska-Bartosik & Torresll2008h would be able to 
account for the stable clock-like X-ray and TeV 7-ray be- 
havior In this scenario, particles move relativistically in the 
pulsar wind, or are accelerated in shocks created by the in- 
teraction of this pulsar wind with the stellar wind. How- 
ever, it should be noted that the H.E.S.S. results do not fa- 
vor this scenario in its standard form, in which the emitter is 
expected to be loc ated between the com pact object and the 
stellar companion ( iBosch-Ramon et all 12008). In this case, 
the absorption is so severe that even cascading cannot help 
to explain disagreement between the measured and predicted 
fluxes around = 0.0. In order to avoid the heavy absorp- 
tion, 7-rays are assumed to be produced at the periphery of 
the binary system inJIakaha shi et alj (l2009l) . similarly as in 
iKhangulvan et al.l(l2008h . 

It is of importance to compare the results with those ob- 
tained for other gamma-ray binaries like PSR B 1259-63 and 
LS I H-61 303. In contrast to the LS 5039 case, the X-ray emis- 
sion from these tw o objects does not seem to show precise 
repea ting features jChernvakova et all 120061 : lUchivama et al.l 
|2009|) . This might be due to the difference of the stellar ob- 
ject type or compact object in the systems. While LS 5039 
host an O star, both PSR B 1259-63 and LS I H-61 303 con- 
tains a Be star, which is characterized with a dense equatorial 
disk. Since a stellar disk plays an important role for X-ray 
emission in the pulsar scenario ( iTavani & Arons|[l997l) . it is 
possible that changes in disk structure avoid short t i me re- 
peating features in X-ray modulation. iRomero et al.l ( l2007h 
actually pointed out that the stellar disk in LS I H-61 303 would 
be likely disrupted around periastron passage by tidal forces 
from the compact object. On the other hand, the emission 
from LS 5039 may be determined only by the location of the 
compact object relative to the O star and thus exhibit precise 
repeating features. 
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